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Abstract 

This paper summarizes preliminary experimental and theoretical research that was 
directed towards the study of quasi-steady-state power flow in a large, un-optimized, 
multi -megawatt coaxial plasma thruster. The report addresses large coaxial thruster 
operation and includes evaluation and interpretation of the experimental results with 
a view to the development of efficient, steady-state megawatt- class magnetoplasma- 
dynamic (MPD) thrusters. 

Experimental studies utilized the Coaxial Thruster Experiment (CTX) facility at the 
Los Alamos National Laboratory. The un-optimized thruster, one meter in length 
with inner and outer tungsten- coated electrode diameters of 0.37 m and 0.56 m 
respectively, was operated over a range of peak power levels (nominally 10 to 40 
MW) in order to ascertain high-power scaling properties. In addition to pure self- 
field operation, an unoptimized applied magnetic field configuration, with adjustable 
flux densities, was used to form a rudimentary, annular magnetic nozzle. Power 
distribution-relevant measurements employed a large complement of diagnostics such 
as spatially resolved C0 2 interferometry, time and space-resolved bolometry, time 
and space resolved electrode IR calorimetry, double-swept Langmuir probes, mag- 
netic probes, CCD optical imaging, and time-of-flight neutral-particle spectroscopy. 
In this initial research, an accurate power balance was precluded due to the lack of 
spatial resolution of power deposition on the electrodes. However, the data, together 
with a calorimetric IR data reduction algorithm, did quantify energy flux deposi- 
tion on electrode surfaces and provided a qualitative picture of global power flow. 
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Estimates of the active electrode surface, which proved to constitute reasonably sub- 
stantial fractions of the total electrode area, indicated that the discharge was well 
distributed over the electrode surfaces. 

In support of the experimental research, a preliminary theoretical model of electrode 
power deposition was developed in conjunction with measured potential profiles. A 
calorimetric data-reduction algorithm was developed to convert the measured tran- 
sient temperature excursions in the electrodes into estimates of peak surface power- 
density deposition from the plasma propellant. Kinetic estimates of the surface 
power-density due to the deposition of charged particle energy acquired by falling 
through the measured anode sheath potential drop, were in semiquantitative agree- 
ment with this algorithm. Fundamental considerations of the plasma production 
process were used to estimate the ion temperature, and by equipartition, the elec- 
tron temperature. The latter was quite comparable to the measured T e , especially 
when allowing for the modest amount of microturbulent heating that may be associ- 
ated with the Lower Hybrid Drift instability under the operating parameters of the 
device. More importantly, the measured longitudinal flow velocity of the propellant 
was in essential agreement with the prediction of self-field nozzle theory, achieving 
Alfvenic flow velocities of ss 10 5 m/s. 

Other highlights of this study indicate that radiative emission is a negligible power 
loss mechanism (< 10%) over the operational power range studied. Also, preliminary 
measurements indicate that the magnetic topology between anode and cathode can 
produce a substantial influence on the electrode sheath potentials, especially on the 
anode fall. This thereby suggests a means for influencing the power deposited on 
the electrodes, with concomitant benefits to thruster efficiency, without relinquishing 
axisymmetry. 
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I. INTRODUCTION 


This report summarizes initial research performed on the Los Alamos Coaxial Thruster 
Experiment (CTX) facility that was directed towards establishing a perspective and 
understanding of quasi-steady-state power utilization in large, high-power, electrody- 
namically driven plasma thrusters. Principal research objectives include developing 
an understanding of thruster efficiency and its dependence on device parameters 
such as spatial dimension, power level, mass-flow rate, and electrode sheath prop- 
erties in the limit of steady-state operation. The ongoing work is motivated by the 
expectation that experience amassed over years of research on coaxial plasma guns 
could be directly applicable to the design of more efficient, steady-state, high-power 
magnetoplasmadvnamic (MPD) thrusters . 1,2 

Coaxial plasma thrusters, like MPD thrusters, utilize a magnetohydrodynamic plasma 
acceleration mechanism to produce thrust . 3,4 In addition to the ionization losses and 
other (nominal) frozen-flow losses, thruster efficiency can be considered as limited by 
characteristic processes occurring in two distinct spatial regimes. One regime, which 
encompasses the part of the annular inter-electrode volume that is dominated by 
the electrode sheath potentials, is characterized by electrode power losses acquired 
by the localized kinetic acceleration of charged particles, and leads to undesirable 
power deposition on the electrodes. This regime has presented a major impediment 
to the development of efficient MPD thrusters . 2 Its study comprises a major focus 
of the work reported here. 

The other regime involves the bulk of the annular volume containing the main flow 
field. As discussed below, very basic considerations regarding the main flow field 
show that it also can have a substantial influence on thruster efficiency, and thus 
its presence and nature is an important constituent of our research. This is a for- 
tiori true when a global power balance is ultimately attempted. Although these 
two spatial regimes are mutually interactive, and may conceivably be separated by 
graduated spatial sub-regimes, this zero-order idealization provides an efficacious 
framework for the study of thruster efficiency. 
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II. the ctx facility and thruster operation 


A photograph and schematic of the CTX coaxial plasma thruster is shown in Figure 1. 
Coaxial guns of this size have been extensively used as a plasma and magnetic helicity 
(linked magnetic flux) sources in magnetic fusion and defense program experiments 
at power levels exceeding 100 MW. 5 ’ 6 Pulsed coaxial plasma guns of similar design, 
albeit with different spatial and operational scales, have been flown in space to 
investigate the dynamic interaction of high-energy plasma with the upper atmosphere 
and the Earth’s geomagnetic field, 7 as well as investigated for advanced propulsion 
applications 8 and materials processing. 9 

For this study, the CTX facility (Figure lc) was configured to operate in a power 
range from 10 MW to 40 MW. 200 A 2 MJ (lOkV) capacitor bank was direct-coupled 
through ignitron switches to the electrodes, where the outer electrode (anode) was set 
at the reference ground potential. A 10 to 40 MW discharge required a bank charge 
voltage of only 1 to 2 kV. At this low charge voltage, the banks accessed only 4% of 
their maximum energy capacity. Direct discharge of the banks into the thruster at 
this low energy, without pulse-forming networks or voltage transformers, produced 
“round-top” current and voltage waveforms of approximately 1 ms in duration as 
shown in Fig. 2. On the magnetohydrodynamic time scale for plasma evolution, 
however, the current and voltage were quasi-static. In this sense, the discharges 
were quasi-steadv-state. 

The CTX thruster consists of two coaxial, stainless steel, tungsten-coated electrodes. 
The inner and outer electrodes are 1 meter long and have diameters of 0.37 m and 
0.56 m, respectively. The working gas, typically hydrogen or deuterium, is injected 
in the annular acceleration region by six fast puff valves that are fired roughly a 
millisecond before the electrodes are energized. The valves are magnetically driven 
in parallel by a current pulse from a small (340 J @ 2 kV) capacitor bank. They 
fully open in roughly 100 fis and empty a 1 cm 3 volume of compressed (up to 10 
atm) gas per valve. On the 1 ms discharge timescale, primary fueling is strongly 
influenced by gas desorbtion off the electrode surfaces. 5 The thruster is fired inside 
a 5 foot (dia) by 15 foot cylindrical vacuum tank. Three 500 1/s turbo- molecular 
pumps and three 25001/s vacuum cryopumps maintain a base pressure of 10 -6 torr. 

200 Operation at power levels up to 900 MW was occasionally accessed for diagnostic calibrations 
and scaling insights. 
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Tank pressure during thruster operation remained less than one millitorr. 

The center electrode (cathode) contains a solenoidal coil that produces an applied B rz 
field (the nozzle field). The coil is driven by an ignitron switched 1.19 mf capacitor 
bank. The nozzle field rises to peak value in 6 ms at which time the discharge is 
initiated. The vacuum field configuration of B rz is approximately constant during 
thruster operation. However, as discussed below, interaction of the flowing plasma 
with the applied field distorts the vacuum field shape and is believed to form a 
rudimentary magnetic nozzle. 5 A calculated estimate of the vacuum B rz field at 
the time of discharge initiation is shown in Figure 3. This unoptimized B rz field 
configuration is uncoventional within MPD thruster research in that the applied 
magnetic field lines are intercepted by both electrodes, with a connection length 
estimated to be not very much larger than the interelectrode spacing (A % 0.1 
m). Thus, from the MPD standpoint, the CTX plasma thruster is unusual with 
regard not only to size and power level, but also with regard to the applied magnetic 
configuration. Future research will focus on optimizing the B rz field configuration 
with respect to optimizing thruster performance. 

Diagnostics, which measure properties within the thruster and the plasma plume, 
included voltage and current sensors, magnetic and electrostatic probes, CCD visible 
imaging, x-ray imaging, infrared electrode calorimetry, multi-chord interferometry, 
visible spectroscopy, neutral particle spectroscopy, and wide bandwidth bolometry. 
This ensemble of diagnostics was chosen to address global power balance including 
radiative, thermal, and particle energy loss channels. Other key facility hardware 
were a computer-controlled safety interlock, bank-charging, and data acquisition 
system. Data was acquired by CAMAC based signal conditioners and digitizers. 

III. QUALITATIVE CHARACTERISTICS EXPECTED OF CTX COAX- 
IAL THRUSTER OPERATION 

Coaxial plasma thrusters have two demonstrable attributes: they produce “high- 
grade” (large magnetic Reynolds number), energetic plasma and they operate over 
a range of spatial scales. As previously discussed, both of these attributes can 
beneficially affect thruster performance. 4 Energetic plasma generation is necessary 
in order to meet the high specific impulse requirements of advanced space missions. 
Furthermore, we believe that efficient steady-state thruster operation requires “high- 
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grade” plasma with good coupling between the plasma and the self and applied 
magnetic fields. Within the context of MHD, these requirements naturally lead to an 
economy of scale for efficient thruster operation with respect to plasma temperature, 
thrust power, and spatial dimension . 4 

Centra] to efficient thruster operation is the utilization of a magnetic nozzle. For 
efficient acceleration of plasmas by magnetic fields, the magnetic pressure gradient 
must push the plasma downstream without slippage. 300 The primary cause of 
such slippage is the presence of resistance to the flow of electrical current in the 
plasma. This electrical resistance constitutes a dissipative property leading to a loss 
of efficiency. To the extent that such resistance can be made small, the plasma 
then approximates the ideal MHD model of a fluid endowed with high electrical 
conductivity. This notion can be quantified in terms of a Magnetic Reynolds Number, 
R = VA/D, where V is the plasma axial flow velocity, A is the width of the annular 
flow channel, and D = rj/p. 0 is the resistive diffusivity of plasma. In the ideal MHD 
model, R is large compared to unity . 4 

The theory of nozzle flow for ideal MHD plasmas has been worked out and provides a 
physical picture somewhat analogous to ordinary gas-dynamic nozzle flow . 10 ' 3 ’ 4 Three 
principal results emerge. First, the acceleration of ideal plasma can be accomplished 
smoothly in a graduated converging-diverging nozzle geometry. Second, the plasma 
reaches a characteristic velocity in the throat of the nozzle. In the absence of thermal 
effects and without an externally applied pervading magnetic field, this velocity is the 
local Alfven speed, Bgj ^fp^p (MKS), where Bg is the azimuthal (“self”) magnetic 
field strength in the throat, p 0 is the magnetic permeability of free space, and p 
is the mass density of plasma in the throat. Third, in the absence of a converging- 
diverging nozzle geometry, there can be no smoothly distributed macroscopic plasma 
acceleration within the context of the axisymmetric, steady-state, ideal MHD model. 
This fundamental result also proves to hold when there is a strong axial magnetic 
field pervading the device . 10 

In the absence of a nozzle, acceleration occurs in nonideal velocity jumps within thin 
dissipative layers, wherein the electrical currents and magnetic pressure gradients 

300 A somewhat similar comment applies to magnetically driven plasma spin-up in the thrusters 
which utilize strong applied magnetic fields. Here the swirl energy is eventually converted to axial 
flow energy in the diverging part of a downstream magnetic nozzle. 
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are concentrated. These layers thus separate uniform regions of differing flow ve- 
locities, and the high current concentrations in these layers can be injurious to the 
electrodes . 301 Moreover, such current sheets are known to be unstable to filamen- 
tation. In contrast, the converging-diverging geometry, in principle, need only rely 
upon a smoothly distributed electrical current to accelerate plasma. 

From the above discussion, we see that MHD plasma accelerators with straight ma- 
terial or magnetic walls must necessarily rely upon localized, nonideal (dissipative) 
properties of plasma in order to achieve plasma acceleration. Hence, MHD accelera- 
tors with straight channels, or with ad-hoc nozzle designs, may be unable to access 
the most efficient coupling of plasma propellant to the driving magnetic fields. For 
example, the Princeton “benchmark thruster” experiments 11 ’ 12 have provided infor- 
mative and detailed data regarding the local concentrations of potential gradients, 
current densities, and heat flux at the lip of a disk-anode that forms a rudimen- 
tary annular nozzle. In this context, Gallimore 12 observed that a severe reduction 
in the nominal anode area (leaving the anode lip exposed) had minimal impact on 
the current- voltage characteristics of the device. This suggests that most of the arc 
discharge is always concentrated at the anode lip. 

A properly designed nozzle geometry can, in principle, produce a smooth, gradual ac- 
celeration of plasma without relying primarily upon the presence of nonideal effects, 
and therefore has access to a more efficient regime of operation. Thus, any investi- 
gation of power flow channels, with a view towards optimization, should be mindful 
of the interrelation between nonideal power losses (including electrode losses) and 
the nature of the main flow field. 

Such considerations are relevant to the applied magnetic field configuration in CTX, 
since the axial flow of electrically conducting plasma imposes a downstream distor- 
tion upon the applied vacuum magnetic field lines. In fact, there is evidence involving 
the “stuffing” and “unstuffing” of the CTX thruster 5 that when the applied magnetic 
field strength is of order or smaller than the nominal strength of the self-field, the 
plasma propellant “blows out” the applied magnetic field. This process forms, at 
least locally, a diverging nozzle configuration just downstream of the muzzle . 5 More- 
over, one may expect that a converging nozzle throat may be formed just upstream 

301 In the less-ideal plasmas, such dissipative acceleration layers would, of course, become thicker 
or more diffuse. 
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of the muzzle. Here, the self-field is “used up” in accelerating the plasma, and hence 
is unable to completely counterbalance the transverse magnetic pressure exerted by 
the applied B rz field. 

These qualitative inferences are supported in CTX by the experimental verification 
of theoretical predictions for self-field-driven flows through converging-diverging an- 
nular nozzles. Measurements of the axial propellant velocity indicate flows at the 
Alfven speed. Because the CTX electrodes are straight, we hypothesize that the 
“nozzle” is provided by the self-consistently distorted applied magnetic field. Direct 
experimental measurement of the nozzle field will be explored in future work. 

To the extent that the present unoptimized CTX thruster produces a magnetic- 
nozzle geometry that is fairly commensurate both with the flow streamlines and 
the equipotential surfaces, one can expect that the electrode effects of acceleration 
and current flow (and concomitant power losses from the plasma) will be gradually 
distributed on the CTX electrodes. This contrasts with thruster configurations where 
the arcs are concentrated in small regions having a constant characteristic location. 12 

IV. ELECTRODE PHENOMENA 

A. Infrared Electrode Calorimetry 

An Inframetrics model 600 IR video camera was used to monitor and quantify surface 
heating on the electrodes. The camera provided absolutely calibrated measurements 
of electrode temperature as discussed in Appendix C. Both end-on and side-on views 
of the cathode and anode were enabled using a combination of salt and germanium 
optics to image the nominal 10 micron infrared emission off the electrode surfaces. 

The Inframetrics camera was operated in either an image mode (17 millisecond 
frame acquisition) or fast line scan (62 microsecond sweep) mode. Both temporal 
and spatial measurements of the heat pulse decay were made. Full co% r erage of the 
electrode surfaces, however, was not possible due to limited optical access. Temper- 
ature excursions from a few degrees Centigrade for 10 MW low power discharges, to 
120 degree increases for 800 MW high power discharges, were observed. These tem- 
perature excursions were converted to peak power deposition estimates by using the 
calorimetric data reduction algorithm discussed in Section IV-B and in Appendix A. 
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In general, shot-to-shot heating nonuniformities were of order a factor of two. This 
variation was attributed to non-reproducible electrode heating due to surface arcing 
or time-dependent, non-azimuthally-symmetric discharge behavior. 

Since the plasma pulse duration and decay of electrode surface heat is short compared 
to the frame acquisition time, the camera, in image mode, gave a distribution of 
temperature increase that was convoluted with the temperature decay during the 
frame acquisition. Thus, in order to resolve the temporal and spatial behavior of 
electrode heating, the camera was operated in the line-scan mode. An example of 
an 4i end-on” infrared image is shown in Figure 4. 

The infrared data can be displayed as either video images or as computer-generated 
plots. An example of digitized line-scan data is shown in Fig. 5a. Here, time runs 
from top to bottom of the frame, and the spatial axis runs left-right. This data was 
obtained by viewing the thruster at a 60 degree angle to the cathode face normal 
from a side-on port. The spatial position of the line scan for this data is shown in 
Fig, 5b. 

In Figure 5a , temperature information during the discharge (CH3.5 ms) is masked 
by noise. However, heat decay after the discharge is clearly visible along the time 
axis. The vertical axis is temperature rise, while the spatial dimension along the 
chosen line-scan is from left to right. Uncertainty in the spatial decay length of the 
temperature is primarily due to the limited viewing access, and to some extent, due 
to the imperfect optical resolution of the IR camera. For the data presented herein, 
these uncertainties precluded absolutely calibrated estimates of total electrode power 
absorbtion. Nevertheless, the IR imaging technique did provide valuable insights into 
the qualitative properties of electrode power deposition as discussed in Section D. 
We note, that this IR electrode calorimetry technique is being refined to allow for 
absolute power deposition measurements in future CTX experiments. 

B. Calorimetric Data Reduction Algorithm 

In this subsection, an algorithm is constructed for converting the observed time- 
dependent surface-temperature increase into an energy flux incident on the local 
electrode surface. This algorithm is based upon two assumptions. First, the radial 
flux of energy out of the plasma is assumed to have the same time-dependence as the 
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electrical input power to the CTX device, which is idealized here as triangular, with 
a linear rise time to peak power, t r , and a linear fall time from peak power, tf. (The 
rise time, t r , shall be taken as the time to reach peak power input to the device; 
the fall time, tj, shall be taken as the time required for the device power to drop 
from its peak value to ten percent thereof.) Second, the tungsten layer of thickness 
1 mm, covering each stainless steel electrode, is taken to be a semi-infinite slab on 
the thermal time scales of interest. This assumption is justified because the thermal 
time constant for a slab of thickness 6 is roughly 0.4 (6 2 /D t h), where D th = (n/pC) 
is the thermal diffusivity of tungsten. 401 This time constant has the value of about 
5.5 ms for 6 = 1 mm. 

For times t small compared to 5.5 ms, we therefore assume that a semi-infinite slab 
of tungsten is irradiated by an energy flux q(t), linearly rising to q max for 0 < t <t r , 
and linearly falling to zero for t r < t < (t r + tf). The problem, then, is to calculate 
the surface temperature rise as a function of time, T(t). Finally, the answer to this 
problem is inverted to express the energy flux q in terms of the observed surface 
temperature rise T. 402 

The heat-flow problem as formulated here can be solved by applying the method of 
Laplace transforms to the heat-diffusion equation. One then finds the temperature 
rise at the plasma-electrode surface to be 

T(r) = qmax /(r,r/) ( ° C) ’ (4J) 

where 


T = t/t r , Tf = tf/t r , (4.2) 

and the function /(r, tj) is given by 

401 Here, k is the thermal conductivity, p is the mass density, and C is the specific heat capacity. 

402 Our approach differs from the previous work of Saber 11 and Gallimore 12 in two respects. First, 
they measured changes in temperature from the back-side of the anode (separated from the plasma by 
the anode wall thickness) whereas we measure temperature excursions on the anode surface directly 
exposed to the plasma. Second, they characterize the energy flux time-dependence as a rectangular 
pulse, consonant with the shape of their flat-topped current pulse. In contrast, the u round-topped” 
CTX discharges, discussed in this report, are best characterized by triangular energy-flux waveforms. 
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f(r,T f ) = u(r) 7 3/2 -u(r-l)(r-l) 3/2 (l+r / 1 )+u(r-l-r / )(r-l-r / ) 3/2 rJ 1 . (4.3) 

Here, u(x) is the unit step-function which vanishes for x < 0 and is unity for x > 0. 

For later reference, the tungsten parameters used in Eq. (4.1) are thermal conduc- 
tivity k, mass density p, and specific heat capacity C, with the standard values 

k — 1.8 x 10 2 (joule/m s °C) 
p = 19 x 10 3 (kg/m 3 ) 

C = 1.3 x 10 2 (joule/kg-°C). (4.4) 

In the present discussion, we provisionally ignore the fact that the sprayed tungsten 
layer probably possesses some porosity, with concomitantly reduced values of mass 
density and thermal conductivity. This will lead us to over-estimate the inferred 
energy flux on the electrode surfaces. 

Examples of the time-dependent surface temperature rise, /(r,r/), are shown in 
Figs. 6 and 7 for t> = 5.5 and 3.0 respectively. A characteristic feature of this time 
dependence is an initially rapid decay transitioning to a slow decay. A similar exper- 
imental result from a nominal 10 MW shot is illustrated in Fig. 8, which is composed 
of a time-sequence of repeated spatial scans, each narrow (60/zs) vertical strip being 
a single spatial scan of the IR camera. Aside from the apparent “spikes” caused by 
the compressed spatial display of each scan, there is some temporal structure during 
the decay, with a 0.25 ms characteristic time scale. This deviation from smooth 
behavior is not understood, but may be associated with the late-time production 
of localized arcs in the gas-starved phase of thruster operation following the main 
pulse. 

In order to utilize Eq. (4.1), to associate the maximum energy flux, q m ar , with the 
maximum observed surface temperature rise, T max , it is necessary to obtain the 
maximum value of /(r, tj). This is shown as the dashed curve in Fig. 9, which is 
also well approximated by the solid straight line, 
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fmax{Tf) % 1 + 0.137/ . 


(4.5) 


Finally, as a rough check, it is also of interest to estimate the time at which the peak 
surface temperature ought to occur according to the present simple model. This 
time (scaled by the rise time t r ) is given by 


r t± ±i£t 

peak " (i + r f y - t) • 


(4.6 - a) 


(By setting r = r peak in /(r,r/), one obtains /max(^/)0 The function r peak is plotted 
in Fig. 10. For r/ > 1, it is well approximated by 


r pca k « 0.49(r/ - 1) + 1.33 . (4.6 - b ) 

Thus, for cases in which the fall time and the rise time are about the same (77 « 1) 
the surface temperature should peak after one third of the decay phase. However, 
for cases in which the fall time greatly exceeds the rise time, 77 ;> 1, the surface 
temperature should peak about halfway through the decay phase. These idealized 
model results appear to be in reasonable qualitative agreement with the observations, 
subject to our earlier remarks that the raw data contain some noise and structure of 
uncertain origin. 

C. Estimates of Plasma Potential and Anode Fall 

The time resolved plasma potential was inferred from electron temperature and float- 
ing potential measurements at three different positions in 40 MW discharges using a 
Langmuir probe. The floating potential during each shot was digitized and recorded 
at the rate of 1 MHz. The probe, originally designed for use as a double, swept 
Langmuir probe, 13 consisted of two 1 mm diameter platinum probe tips separated 
by 3 mm and extending 1 mm beyond a boron nitride insulating mount. The plasma 
potential was obtained from the floating potential using the standard relation: 14 

= tv* + i) In (g|) , (4.7) 
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where kT e was nominally measured at 15 eV in deuterium. 404 

Fig. 11 shows measurements of plasma potential made at three probe positions: 
(1) inside the muzzle, several millimeters from the anode; (2) inside the muzzle, 
5 cm from the anode; and (3) 40 cm downstream from the exit plane. The data 
exhibit high frequency oscillations as well as long timescale variations. Shot to shot 
variations were less than 50 volts. 

Power spectra of the voltage, plasma potential within the muzzle, and interferometer 
signal (center chord, 40 cm downstream of the exit) are shown in Fig. 12. These 
power spectra were ensemble averaged over 8 shots and have a nominal 250 kHz 
bandwidth. The oscillations in the discharge voltage and plasma potential are broad 
band, with a broad peak near 100 kHz. This spectral feature is too low for fundamen- 
tal plasma oscillations. However, since this frequency corresponds to the reciprocal 
axial magnetosonic transit time, it is probably related to macroscopic plasma flow 
dynamics. 

With respect to the slow variation of the plasma potential, the data exhibit two dis- 
tinct behaviors. The plasma potential is substantially positive relative to the anode 
for the first 200 /zs of the discharge, indicating a significantly reversed anode fall 
around the time of peak power. Following a short transition period, the potential 
becomes negative for times greater than 250 fis. The probe floating potential, mea- 
sured at probe positions inside the muzzle, was averaged over time for these two 
periods and these results were averaged over an ensemble of shots taken under the 
same conditions. These results, along with the anode and cathode potentials, yield 
the potential profile in the muzzle (Fig. 13) for two different periods of the discharge. 
The plasma potential several millimeters from the anode (Fig. 13) is positive relative 
to the anode (confining plasma electrons) during the first 200 fis of the discharge 
and negative for the remainder of the discharge. 

The long timescale variation in the plasma potential, measured in the muzzle (Fig 11), 
can be qualitatively explained by investigating the evolution of the magnetic field 
structure during a discharge. Calculation of the vacuum magnetic field in the thruster 

404 The collisionless, thin-sheath approximation with unmagnetized ions is adequate for analyzing 
CTX conditions . 15 A clean probe surface was ensured by arc-discharge cleaning within the vacuum 
vessel. 
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at the time of a discharge is shown in Fig. 3. At discharge initiation, field lines from 
the coil within the cathode provide a magnetic connection from cathode to anode. 
Electrons are only marginally collisional (for classical Coulomb collisions) on the 
spatial scale of the CTX thruster. Hence, we expect the arc voltage to acceler- 
ate current carrying electrons to moderately high energies with respect to the bulk 
plasma electrons; forming a high-energy tail on the electron distribution function. 

405 

Early in the discharge, the plasma potential several millimeters from the anode is 
~225 V positive, confining the energetic electrons. As the discharge evolves, the field 
lines, which connect the cathode to the anode at early times, are stretched by the 
plasma and eventually magnetically connect the cathode to the tank wall. During 
this latter phase of the discharge, the magnetic connection between cathode and 
anode is reduced, and results in an “electron-attracting” anode sheath. 

This behavior is evident in Fig. 14, which shows both the fraction of the discharge 
current flowing to the tank wall and the plasma potential in the muzzle. During the 
early portion of the discharge (< 200 /rs), the plasma potential several millimeters 
from the anode is positive (confining electrons) and the fraction of the discharge 
current flowing to the tank wall is «30%. For times greater than 250 ^xs, the fraction 
of current flowing to the tank wall is «80% and the plasma potential changes sign, 
becoming negative (normal MPD anode fall). 

The behavior of the plasma potential as a function of applied magnetic field structure 
shows that the presence of such an applied field can exert a very strong influence 
on the anode fall voltage. In fact, the anode fall can actually be reversed by a 
substantial amount. This effect is reminiscent of the use of external magnetic cusp 
fields to control the effective anode area and anode fall in ion sources. 16 However, in 
contrast to the ion source results, the magnetic configuration in CTX affected the 
anode fall voltage apparently without a marked loss of axisymmetry. The implication 
is that there may exist an optimal axisymmetric applied magnetic field configuration 
that minimizes the electrode fall voltages and hence minimizes the electrode power 
deposition. 

<05 tVe note that a complete treatment of this physical process must include the effects of micro- 
turbulent beam-plasma interactions. 
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D. Power Deposition Estimates 

Inverting Eq. (4.1), we have at the time when T(°C) peaks out, 


Qmax — 


j KpC r moI (°c ’) 

V ^ r fmax 


(4.8) 


With the standard parameters for tungsten from Eq. (4.4), this expression becomes 


Qma 


joules 


s-m^ 


2.8 x 10 6 T max {°C ) 


V%.(100/is) fmax (■ T f ) 


or also, 


Qma 


*MW1 = 2.8 T max (°C) 

. m^ . yj \ £ r ( 100/is) fmax(l~f) 


(4.9) 


Here, Z r (100/is) is the rise time of the power input to the plasma thruster expressed 
in units of 100 /is. In addition, we recall that f ma x( T j) is approximately a linear 
function of the ratio of the decay time to the rise time, given by Eq. (4.5) with 
T f — tfjt r . In principle, q m ax and T max in Eq. (4.9) are to be interpreted as peak 
values in time for each position on the electrode surfaces. In practice, we utilize 
measurements just after the shot, at around 0.5 ms, to avoid electrical interference. 


The procedure discussed above, hereafter called the first procedure, assumes that the 
surface temperature excursion peaks at 0.5 ms. The results from our data ensemble 
based on this assumption are presented in Tables I and II. In reality, this assump- 
tion probably overconstrains the data ensemble. Thus, we also utilize an alternate 
procedure, hereafter called the second procedure, that allows for the temperature 
excursion to peak at an earlier time in the discharge. The second procedure uses Eq. 
(4.1) with T measured at 0.5 ms and /(r, tj) evaluated at 0.5 ms without invoking 
any systematic offsets or time delays. Results from the second procedure are pre- 
sented in Tables IA and IIA. We believe that these two procedures will bracket the 
energy flux estimates, with the first and second procedures yielding under and over 
estimates, respectively. 


As an example of the first procedure, we consider shot 20383, a 10 MW class dis- 
charge. From Fig. 15, we see that the input power rises to about 15 MW in roughly 
100 /is. Moreover, the main decay of this power is essentially complete (down to 
10% of peak power) by 350 /is. Thus t r % 100 /is, tj « 250 /is, and tj zz 2.5. From 
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Eq. (4.5), one then would estimate that f max wl.33. Fig. 5b, shows a schematic 
drawing of the line scan of the IR camera view for this shot. The scan passes just 
under the cathode, allowing an axial scan of the inner-anode far wall that is estimated 
to cover about 30 cm of length. Beyond this point the view upstream is occluded 
by the outer-anode near wall. In Figure 5a, a space and time-resolved record of 
the temperature rise on the anode is shown. Traversing from left to right at late 
times (6 ms), the lip of the far anode wall is clearly identified by the temperature 
increment at scan position 0.25, and occlusion by the near anode wall is likewise 
identified by the temperature decrement at scan position 0.65. The temporal and 
spatial temperature record is believed to be reliable for times exceeding 0.5 ms. 

From Fig. 5a, we estimate the peak (vs position) temperature rise at about t = 
0.5 ms to be ss 7.7 °C at the lip of the anode, exhibiting a spatial decrease upstream 
along the anode. (Here, 0.47 on the vertical axis is calibrated to 10°C.) The knee 
of the curve of T(°C) vs scan position (due to occlusion by the near anode wall 
in association with the limited resolution of the optics) occurs at about 0.65. The 
temperature rise at the knee is estimated to be about 5.4°C. Using Eq. (4.9), we 
thus estimate the peak (in time) energy fluxes to be 

MW 

Qmax (anode lip) « 15 — =- 

m/ 

MW 

Qmax (“knee” at scan pos. 0.65) « 11 — -x- (4.10) 

nr 

The axial distance from the anode lip to the position of the temperature “knee” 
is estimated to be about 25 cm, from which one can linearly extrapolate to a zero 
energy flux at L ss 85 cm. Using the peak q max value in Eq. (4.10), one can estimate 
the power deposition on the anode (Anode area = 1.76 m 2 ) to be given by (0.5 x 15 
x 0.85 x 1.76) MW, or 

P,w e (20383) w 11MW . (/, = 85cm.) (4.11) 

If a like amount of power is deposited on the cathode (for which there is no infor- 
mation regarding this particular shot), then the power to the electrodes exceeds, by 
roughly 50%, the input power of 15 MW. This overestimate in total electrode power 
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absorption demonstrates the need for better spatial resolution of the electrode energy 
flux distribution. 

In Tables I and II, we summarize the IR imaging parameters for an ensemble of 10 
MW class and 40 MW class discharges. For each respective class, all experimentally 
controlled parameters (i.e. bank charge voltages, gas valve plenum pressure and tim- 
ing, B rz field values and timing) were fixed. The variations in measured parameters, 
such as peak discharge power and g max , reflect statistical variations. Ensemble av- 
eraging over the discharges yielded, for the 10 MW class ensemble, values of average 
peak discharge power and q ma x °f H i 2 MW and 12 ± 3 MW/m 2 , respectively. 
For the 40 MW ensemble, averaged peak discharge power and q ma x were 48 ± 5 MW 
and 23 i 5 MW/m 2 , respectively. Clearly, for both 10 and 40 MW class shots, the 
entire electrode surface is not absorbing the peak energy flux. Damage patterns on 
the electrode surfaces suggest that roughly 25% of the total electrode area of 2.9 m 2 
is involved in high energy flux absorption. 


E. Interpretation 

The standard model for anode energy deposition is based upon the use of fluid- type 
energy-conservation equations for electrons and ions. 


d_ 

dt 


(flow-energy density -I- thermal-energy density) 


+ V • (enthalpy flux + viscous power flux + heat flux) 


= rate of work per unit volume performed by the electric field. (4-14) 

Here, the enthalpy flux is the sum of the flow and thermal energy density fluxes, 
plus the rate of pressure- work per unit volume performed by an element of fluid on 
the surrounding fluid; and the viscous power flux comes from the rate of work done 
by viscous forces in the velocity-shear layer. (Collision terms that ultimately cancel 
upon addition of the electron and ion equations have been omitted.) Orienting the 
x axis (radial direction) into the anode, one assumes a steady-state, 1-D model, and 
integrates across the anode fall from the “plasma edge” (x = x p ) up to the surface of 
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the anode (x = x Q ). Adding the resulting electron and ion equations, but neglecting 
electron inertial and viscous terms, one finds the power density on the anode, P Anode 
to be expressed as (with 7 = 5/3) 


P Anode — + § (Ti + T e )) + <}ix + 9ex + ion viscous power flux] 

j x=ip 


5 

2 


Jx 


T 

— + J X $ 

e 


p 


+ work function term + plasma radiation term , 


(4.15) 


where the additional effects of the anode work function and radiation from the plasma 
have been explicitly inserted. Here, m is the ion number density, V t is the ion 
fluid velocity, £( e j is the ion (electron) heat flux vector, J x (< 0) is the uniform 
radial current density, Tq e ) is the ion (electron) temperature at x = x p , and 
is the plasma potential relative to the anode at the nearby plasma edge, x = x 
Thus, when < 0, electrons are accelerated into the anode. Abbreviated forms 
of Eq. (4.15) are sometimes employed to estimate electron power deposition on the 
anode in MPD thrusters. 11,12 

We maintain that the fluid description epitomized by Eq. (4.15) is inappropriate 
for describing electrode power absorbtion in CTX. The Coulomb-collision mean free 
path in CTX at the nominal 40 MW level is many centimeters. Moreover, the 
deuterium ion gyro-radius is at least of order 1 cm and the ions are highly magnetized 
( u, « ~ lO^ti). We conclude that for electrode or sheath phenomena, the plasma 
propellant is highly kinetic in nature. Therefore, it is unrealistic and unreliable to 
employ a fluid model to describe physical processes that are occurring over kinetic 
distances. Were such an attempt to be made, one would immediately encounter 
the problem that the generally important thermal conduction and viscous terms 
in Eq. (4.15) become large, ill-defined, and not usable. Moreover, their use would 
require some prior knowledge of the thermal and velocity gradients in the vicinity of 
the near- anode boundary layer. 

Instead, we employ a simple and fundamentally kinetic estimate of charged particle 
power deposition onto the anode for 40 MW class shots. We note that during peak 
power, these shots actually exhibit a reversed anode fall, wherein the plasma is 
positive relative to the anode by as much as 220 V (time-averaged) within several 
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millimeters of the anode surface. Thus, it is expected that plasma ions primarily are 
responsible for the observed power deposition during the peak power phase of these 
shots. The bulk-plasma electrons, with temperatures measured between 10 eV and 
20 eV, do not play a significant role in anode heating. 

The ion power deposition onto the anode in 40-MW shots can be estimated by 
calculating the ion thermal flux from the edge of the plasma towards the anode, 
and taking into account that the ions will acquire about 200 eV in falling through 
the reversed anode sheath. The exact ion temperature will little affect estimates 
of this ion thermal flux, and moreover, will prove to be only a minor addition to 
the 200 eV potential fall. The axial ion kinetic energy (~100 eV) is ignored in 
constructing this estimate because the radial profile of axial velocity is not well 
known in this device, and therefore it is possible that a boundary layer of reduced 
axial velocity exists adjacent to the anode. This approach will, therefore, provide a 
modest underestimate of the ion power deposition. 

Although the axial flow energy in the main flow field will not be directly convected 
to the wall in a no-slip viscous boundary layer, it does supply energy for viscous 
heating within that layer. This extra heat will be partly conducted and convected 
downstream. Accurate estimates of these various effects in the highly kinetic regime 
would require the solution of kinetic equations for the non- Maxwellian deformations 
of the ion distribution function. We note that these neglected non-Maxwellian effects 
are driven by ion energies that do not dominate the reversed anode fall voltage, 
and probably contribute to only moderate increases in the estimates of ion power 
deposition onto the anode. 

Ion temperature has not been directly measured within the thruster. Instead, it 
is estimated here from fundamental considerations of the plasma- production and 
entrainment process (see Appendix B). The T,-, so estimated, proves to yield semi- 
quantitative agreement with certain measurements of transverse energy in the plume. 
The fundamental observation here is that when a given mass at rest is abruptly leg- 
islated into a moving system of a specified velocity, simultaneous conservation of 
momentum and energy forces the production of heat, such that the flow energy and 
the thermal energy are equalized. A simple example of this phenomenon in gas dy- 
namics is the state of a shocked gas behind a piston-driven strong shock propagating 
into a gas at rest. One finds there that the post-shock thermal energy density equals 
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the flow energy density. Similarly, ions created from neutrals at rest or moving slowly 
relative to the propellant in the coaxial plasma thruster suddenly find themselves 
in crossed E and B fields, forced to move at the E/B drift velocity. Irrespective 
of collisionality, it can then be shown that the ion thermal energy and drift kinetic 
energy must be equal. 4 (This effect is sometimes referred to as an “ion-slip” effect.) 

The “stuffing vs unstuffing” property of CTX thruster operation leads us to infer 
that the relevant drift velocity seen by newly created ions is ( E r /Bg ), since the 
applied magnetic field must be pushed aside by the self-field magneto-plasma, and 
because independent measurements yield the result that B 2 < Bj. For 40 MW 
shots, typical values are near E r 4 x 10 3 V/m and Be ~ 0.08 Tesla in mid-channel, 
leading to initial drift velocities (created in the upstream region) of about 0.5 x 10 5 
m/s. For deuterium, this value corresponds to drift energies, and hence thermal 
energies, of about 25 eV. (Ions created near the nozzle throat at the muzzle would 
have drift velocities, hence thermal velocities, near the Alfven speed, but would be 
swept downstream away from the muzzle.) Equating the ion thermal velocity to the 
drift velocity of newly-created ions, and using a representative density measurement, 
n as 1.0 X 10 20 m -3 for 40 MW shots, we arrive at an estimate of the ion thermal 
flux, T t ^j. Using the canonical expression for T t /u, 850 

_ 1 1 „ 9 r particles 

r thx ~ -nvthi = - x 10 25 - — j . (4.16) 

4 8 m -s v ’ 

Since each ion has acquired about 200 eV in falling through a gyro-radius spanning 
the reversed anode sheath, we arrive at an estimate of the ion power density at the 
anode, namely, 


P. Anode ~ r t/l , x 200 eV = 40 MW/m 2 . (4.17) 

This result is in semi-quantitative agreement with the 24 to 31 MW/m 2 average 
energy flux shown in Tables I and la. 

In summary, an algorithm has been developed to convert the observed transient 
surface heating of the electrodes into estimates of incident peak surface power density, 

SS0 The expression nv t h,/(2^/P) would be slightly more accurate for a half-Maxwellian distribution 
function. 
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and these estimates are in semi quantitative agreement with a simple model of ion 
thermal flux into the anode during the peak power phase with a reversed anode fall 
voltage. 

V. GLOBAL DISTRIBUTIONS OF POWER 
A. Radiation 

A new-style windowless XUV photodiode, absolutely calibrated at the University of 
Washington, was used with a scanned collimator to measure the absolute radiation 
losses from the end of the CTX plasma thruster. The bolometer consisted of an 
XUV100C silicon photodiode detector, current-to-voltage preamp, fiber-optic link, 
vacuum housing, and pinhole aperture and collimator. These detectors have a highly 
stable internal quantum efficiency and flat spectral response up to 6 keV energy pho- 
tons. 1 7 The XUV photodiode sensitivity is 0.275 Amps/Watt (radiation incident). 
The bolometer has been cross calibrated against a relatively traditional (but more 
insensitive and bulky) gold foil calorimeter from Los Alamos, 18 which was itself cal- 
ibrated using a KrF pulsed laser at 248 nanometer wavelength. 19 The advantage of 
the XUV photodiode lies in its simplicity, greater sensitivity, and speed of response. 
Its principle disadvantage is the uncharacterized response to possible particle fluxes. 

The bolometer viewed the plasma through a 100 micron diameter pinhole aperture, 
which was also limited by a 12 cm long, 0.79 cm diameter threaded tube to define 
the field of view. By taking into account the detector solid angle, correcting for 
the fraction of the plasma viewed, and assuming an isotropic radiation source, we 
calculate that only a tiny fraction of the plasma radiation (2.3 • 10~ n )) hit the 
detector. In calibration, 1 /zAmp out of the detector corresponds to 700 kW of 
total plasma radiation, with estimated 20% absolute error bars. Noise levels were 
at 100 nanoamps, or lower, and the instrumental response time (with electronics 
package) was better than 1 /zs, but limited to 4 /zs to avoid aliasing at the digitizer. 

By scanning the collimated bolometer across an end-on view of the thruster, on 
a shot-to-shot basis, we observed that the radiation was largely confined spatially 
to the anode-cathode gap (Figure 16). If one assumes isotropic emission of the 
radiation, and cylindrical symmetry of the emission pattern, w T e then infer an absolute 
radiative power loss of 3-6% ±1% of the peak I-V input power to the thruster, 
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when operating over the 10-40 Megawatt range. In the case of extreme high power 
operation (« 900 MW), the radiation fraction was observed to increase, although 
only to 10% of the input power. In summary, absolute calibrated bolometry gives 
an upper bound on the radiation losses, which range from 3-6% for 10-40 MW input 
power. We therefore conclude that radiation loss is a small constituent of global 
power flow for the modes of operation studied. 

B. Estimates of Thruster Performance 

A double-swept Langmuir probe (DSP), a CO 2 multichord interferometer (MCI), 
and a time-of- flight (TOF) neutral particle spectrometer were used to ascertain the 
intrinsic properties of magnetized thruster-generated plasma flow. Included were 
measurements of electron temperature and ion density in the muzzle and plasma 
plume, and direct measurements of axial (longitudinal) neutral particle flow in the 
plume. In this section, we focus on thruster performance during nominal 40 MW 
operation. 4 

From spatially-resolved multi-chord CO 2 interferometry we infer that the ejected 
plasma plume is of order the thruster diameter at an axial position 40 cm downstream 
of the muzzle. This result is roughly independent of input power and depends mainly 
on the ratio of the applied (B r2 ) to self (B$) magnetic field for a given discharge. 
Self-field operation generally results in a density distribution that is peaked on axis. 
The density distribution broadens as a function of the applied to self-field ratio until 
the thruster reaches a limit where plasma ejection terminates. Qualitatively, this 
“stuffing limit” is reached when the tension along the applied B rz field effectively 
exceeds the pressure in the self B$ field and prevents the net ejection of plasma. This 
effect is demonstrable over the entire power range of thruster operation. 4 

During nominal 40 MW operation, probe measurements in the plume 40 cm down- 
stream of the muzzle at the interferometer plane suggest a mean ion density of 
« 3 ± 1 • 10 19 m~ 3 . Reasonable agreement exists between the DSP and MCI density 
estimates when the plume profile is of order the thruster diameter. DSP inferred 
electron temperature in the plume is between 10 eV and 15 eV. 

The DSP was utilized to compare plasma parameters in the plume and muzzle at a 
fixed radial position of r « 0.25 m (referenced to the thruster axis of symmetry). This 
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radial position placed the probe roughly half-way between the anode and cathode 
when inside the muzzle. Both temperature and density measurements in the muzzle 
were complicated by plasma noise and probe arcing. Nevertheless, a limited set 
of data provided estimates (within a factor of 2) of these parameters. The mean 
density at the muzzle was 3 to 4 times the downstream density. Again, for nominal 
40 MW operation, the density at the muzzle was estimated at « 1.0 ± 0.5 * 10 2 ° 
m -3 . The electron temperature in the muzzle was « 10 to 20 eV, similar to the 
value downstream. 

Measurements of plasma density in the muzzle, and concomitant projections of den- 
sity in the breech (from a Bernoulli flow model) allow preliminary comparisons to 
be made with the predicted MHD performance. 

Plasma Flow From our previously discussed model of plasma flow in a converging- 
diverging magnetic nozzle, the plasma flow should be Alfven like (V e ~ V^CAb)- 4 
Measurements of the axial neutral flow were made with the time-of-flight spectrom- 
eter. An example of a TOF arrival time distribution is shown in Fig. 17. The 
difference between the first signal spike at 135 /xs (the timing fiducial) and the sec- 
ond large signal spike at 153 jus (arrival of the neutral particle flux) corresponds to 
an axial flow velocity of « 1.0 ± 0.3 • 10 5 m/s. Furthermore, the narrow shape of the 
spike is indicative of an energy distribution where the directed flow energy is larger 
than the longitudinal thermal energy. 4 

In comparison, the estimated Alfven velocity, based on density measurements in the 
muzzle, was Cao ~ 0.9 ± 0.2 • 10 5 m/s. The MHD model predicts a flow velocity for 
fully accelerated ions of V e = y/ZC A o ~ 1-5 ±0.4 -10 s m/s, 4 and is in reasonably good 
agreement with the neutral particle velocity data. In contrast, the E x B drift ve- 
locity for this discharge has an upper limit in the breech of Veb < E r /Be ~ 0.5 • 10 3 
m/s. In general, the axial neutral flow velocity is observed to exceed the E x B drift 
velocity in the breech by a factor of 2 to 3 early in the discharge when the thrustei 
is far from the stuffing limit. The basic agreement between the measured flow ve- 
locity and the MHD model prediction suggests that converging-diverging nozzle flow 
may be an important constituent of gross plasma acceleration in the present CTX 
thruster. 

Plasma Quality and Electrical Effort The measured values of plasma density, 
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temperature, and flow velocity allow an assessment of “plasma quality”. Assuming 
classical resistivity, these values project a magnetic Reynolds number of approxi- 
mately 1300, indicating good coupling between the plasma propellant and the mag- 
netic fields. 

Another important operational parameter is the number of elementary electrical 
charges that must be transferred between the coaxial electrodes for each ion-electron 
pair of propellant that is ejected downstream (for thrust). This quantity, which we 
term the Electrical Effort E, may be expressed as a dimensionless number, 

~ = (mi/e)(I/M) . ( 5 . 1 ) 

where m t - is the mass of a propellant ion, e is the elementary electron charge, M is the 
mass flow, and I is the discharge current. A large value for the electrical effort means 
that large discharge currents are required to process a given mass throughput, with 
concomitant wear and tear on the electrodes. Thus, the reciprocal of the electrical 
effort can be thought of as a figure of merit for any electrodynamic thruster. 

If the measured neutral particle flow velocity accurately indicates the plasma flow 
velocity at the muzzle, the estimated mass flow, based on the density and flow 
velocity measurements, is M = 27rp m V m r m A « 4.5 gm/s. Here, p m and V m are the 
mean plasma density and flow velocity in the muzzle, respectively, r m is the mid- 
channel radius, and A is the channel width. This gives an electrical effort of E % 0.5. 
These values of Rm and E are indicative of fairly high-grade plasma and consistent 
with the observed ideal MHD-like thruster performance. 

Thrust Power and Efficiency With all due regard to relevant uncertainties 
in the experimental measurements and technique, it is interesting to note that the 
combination of the measured axial flow velocity and the estimated mass flow rate 
yields a peak thrust power estimate of 23 MW. The corresponding peak thrust for 
these conditions is 500 Nt. This value is close to that calculated from the MHD 
model which predicts an upper thrust limit of 580 Nt. 4 

Thruster efficiency, 6, can be defined in terms of the electrical effort, E, as 


where e z is the mean axial energy of the propellant (in volts) and Vj is the electri- 
cal discharge voltage that appears across the electrodes of the thruster. From this 
expression, we conclude that the electrical effort must always exceed the ratio of 
the ejected energy per particle (in volts) to the discharge voltage. For the 40 MW 
class discharges, our 0.5 estimate for E, coupled with a measured axial flow energy 
of « 100 eV (for deuterium) and nominal discharge voltages of « 400 V, projects a 
thrust production efficiency of « 50 %. This value lies within the range of efficiencies 
reported for pulsed, high power coaxial plasma thruster operation in xenon. 8 

C. Thermal Plume Losses 

Plasma Heating As previously discussed, ion entrainment in a flow field coupled 
with ion- electron equipartition results in plasma heating. The observed 10 eV to 
20 eV electron temperature in the muzzle is understandable from this perspective. 
Over the 10 MW to 40 MW discharge conditions, ion temperature is predicted to 
reside in a 25 eV to 45 eV range as compared to axial flow energies of 100 eV to 200 
eV (see Appendix B). However, for energies exceeding 20 eV, the electrons effectively 
decouple from the ions as their classical equipartition time exceeds their axial transit 
time. Hence, the observations of electron temperature in the 10 eV to 20 eV range 
are consistent with heating to their equipartition limit by the higher temperature 
ions. They are also consistent with the estimated ineffectiveness of classical ohmic 
heating on the timescale of an electron axial transit. 4 It is possible that a moderate 
contribution to the electron temperature is provided by microturbulence associated 
with the Lower Hybrid Drift and related instabilities; and by heating due to collective 
interactions with primary electrons from the cathode. These issues will be discussed 
in future work. 

VI. SUMMARY 

In summary, various measurements and associated theoretical interpretations have 
been implemented in a preliminary investigation of power-flow in the CTX plasma 
thruster. The results are, so far, constrained by the limited number of measurements 
taken, the uncertainty inherent in certain measurements (e.g., the active electrode 
area), the lack of instantaneous mass flow control, the unoptimized electrical circuitry 
driving the thruster, and the unoptimized nature of the applied magnetic field config- 
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uration. These constraints will be addressed in future research by ongoing upgrades 
in the electrical circuity and mass flow control that will drive flat-topped, long-pulse 
(« 10 ms) discharges at 1 to 40 MW. 

In our preliminary investigation, good agreement was found between the measured 
properties of CTX thruster operation and the predictions of a model based on ideal 
MHD. This agreement suggests that ideal MHD processes play an important role 
in coaxial plasma thruster dynamics. Within this context, CTX performance with 
respect to propellant flow velocity, thrust production, and operational impedance are 
reasonably well understood. Interpretation of the preliminary data also reinforces our 
belief that optimized magnetic nozzle design can be beneficial to thruster efficiency. 
Related considerations carried on in parallel to this work have lead to the conceptual 
design of an electrically adjustable, axisymmetric, annular magnetic nozzle that is 
intended to address the simultaneous control of plasma acceleration, electrode sheath 
potentials, and plasma detachment. 

It is likely that the perceptions and concepts presented in this paper have efficacy in 
the development of efficient, megawatt-class MPD thrusters. Specifically, designing 
MPD thrusters (in power and spatial scale) to access efficient plasma-magnetic field 
coupling and utilizing an optimized B rz field configuration to maximize efficiency 
and minimize anode fall could prove advantageous. (We note that these perceptions 
are consonant with previous quasi-one-dimensional analysis of flow channel geometry 
effects on MPD thruster performance. 20 ) Future coaxial thruster research will focus 
on these issues in the context of scaling multi-megawatt thruster performance to the 
smaller spatial scales and operational power levels required for MPD devices. 
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Appendix A 

Derivation of the Algorithm for Converting 
Surface Temperature into Incident Energy Flux. 

The energy conservation equation that governs heat diffusion in a material 
free of internal heat sources (and sinks) is 

dr 

P C -fa =-V-?, (A-l) 

where T is the local temperature, p is the local mass density, C is the heat 
capacity per unit mass (specific heat capacity), and q is the local internal heat 
flux. 

This local heat flux, q, is related to the local temperature gradient, VT. by 


q - -kVT , (A-2) 

where k is the local thermal conductivity of the material. 

Combining (A-l) and (A-2) assuming uniform k, one obtains the equation 
for temperature diffusion in a one-dimensional slab model, namely, 


dT d 2 T 
dt ~ Dth dx 2 ’ 

wherein the thermal diffusivity of the material, D t h , is defined by 


(A-3) 




(A-l) 


We shall henceforth assume that p (as well as k) is uniform, so that D t h is 
uniform. Concomitantly, we shall assume that k, p. and C are all independent 
of the (nonuniform) temperature, which is a good approximation over the lim- 
ited range of temperature excursions that were observed on the CTX electrodes 
(10°C-20°C, for the 10-40MW operating range). Moreover, we note that the 
symbol, T, in Eq. (A-3), obviously can be interpreted as the temperature ex- 
cursion away from the initial, uniform temperature; or more generally, as the 
temperature difference from any convenient reference temperature. 

We shall now characterize the thickness of the tungsten electrode coating 
(s: 1 mm) within the context of a heat diffusion problem for which the surface 
irradiation time is very short. The intrinsic thermal time constant for a slab of 
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thickness 6 having thermal diffusivity D t h, is of order 6 2 /D t h , since <5 and D th 
are the only characteristic quantities in the diffusion problem when the external 
conditions are changed very rapidly compared to 6 2 /D t h. Additional detailed 
considerations sketched here show that the effective thermal time constant is 
more accurately represented by about 0.4<$ 2 /Z? t h- For example, consider the 
situation in which a slab is initially at some given uniform temperature, but 
for which the temperature of one side (the "near” side) is suddenly lowered 
(by external forcing) to a new value and then held there. (If the near-side 
temperature were to be suddenly raised instead of lowered, the problem could 
be trivially reformulated so as to be identical to the present case.) Assuming that 
the far side of the slab is insulated, how long will it take for the temperature of 
the far side to approach the new temperature of the ’‘near'’ side? This problem 
may be solved by expanding the function of position, x, and time, t, namely 
T(x,t), in a spatial cosine series with time-dependent coefficients, each term of 
which satisfies both boundary conditions. Using Eq. (A-3). one finds that the 
slowest-decaying mode is the cosine term having only one-quarter wavelenth in 
the slab, for which the decay time is (4/7r 2 )(d 2 /D th ). This global mode makes 
the major contribution to the decay of the temperature at the far side. (The 
higher modes, which are of significantly smaller magnitude, are more relevant to 
the decay of “sharpness” of the square corner of the temperature profile at the 
near side.) This elementary viewpoint of the thermal time constant of a slab is 
borne out by numerical solutions. 

For example, consider the problem more relevant to our situation, a slab 
initially at some uniform temperature, suddenly irradiated at time t = 0 by a 
constant heat flux at the near side (q(t) = 0 for t < 0, q{t) = q 0 for t > 0). 
Detailed numerical solutions of the time-dependent temperature profile for this 
particular problem, in which the far side is insulated, can be found in the book 
Unified Analysis and Solutions of Heat and Mass Diffusion by M. D. Mikhailov 
and M. N. Ozi§ik, John Wiley and Sons, 1984. Fig. 7.4-a, page 302. There, 
normalized time r is expressed in units of (6 2 /D th ). One sees clearly from 
this figure that, for times r<0.4, the temperature of the far side remains very 
small compared to the temperature of the irradiated side; whereas for increasing 
r, r > 0.4. the continuously increasing temperatures of the two sides become 
comparable. (Detailed numerical solutions to the preceding problem of a sud- 


30 



denly lowered surface temperature are presented in Fig. 7.1-a on page 289 of the 
above reference. One again finds that r = 0.4 constitutes a reasonable charac- 
teristic time beyond which the "far” side is influenced by the new temperature 
of the “near” side.) Thus 0AS 2 /D t } l constitutes a reasonable time-demarcation 
for the effective thermal time constant for heat diffusion in a slab. 

The conclusion that can be drawn from the above discussion is that, for 
surface irradiation times (and observation times) short compared with QAS 2 /D t h, 
one can effectively regard the slab as being very thick (semi-infinite). Since 
the thermal time constant for a 1-mm-thick tungsten layer is about 5 ms, and 
since the CTX pulse duration is £1 ms, and since our observation times will 
be focused around 0.5-1 ms from “turn-on” time, we shall justifiably invoke the 
simplification that the tungsten layer is indeed semi-infinite within the context of 
highly transient heat diffusion induced by a short pulse of one-sided irradiation. 

Using the above simplification, we shall now derive the solution for the time- 
dependent surface temperature rise induced by a time-dependent energy flux 
incident on the surface. The method used is the application of the Laplace 
transform to the temperature diffusion equation, with the use of appropriate 
boundary conditions. 

Define a normalized time r by r = D^t (with the dimension of length- 
squared). Then Eq. (A-3) becomes 

dT d 2 T . 

dr ~ dx 2 ' ^ A ' 5 ) 

This diffusion equation is to be solved for the surface temperature excursion (at 
x = 0), subject to the initial condition that 

T{x, 0) = 0 , (A-6) 

(in which T — 0 consitutes the ambient reference temperature level); and the 
boundary condition at the “near” side, 

fdT\ 1 , s 

(& L„ = -p {r) ■ (A - 7) 

wherein the time-dependent energy flux, q(T), incident on the near surface at 
x — 0, will be specified later below. The boundary condition at the very distant 
boundary (the "far" side) will be taken as 
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T( oo. r) = 0 , 


(A-8) 


in accordance with the transient approximation discussed above. 

Recalling that the Laplace transform is defined by 

0(x,p) = C p [r(x, r)] = f T(x, T')e~ VT ' dr' . 

J o 

wherein Re(p ) > 0. one takes the transform of Eq. (A-o) (noting the initial 
condition Eq. (A-6)) and finds 


d? = pe ’ < A -«) 

with the general solution 

0(x. p) = A p e x ^ + B p e~ x ^ . ( A-10) 

Without loss of generality, we shall use the convention that when Re{p) > 0, 
Re(Vp) > 0. At far distant positions, the temperature excursion T(x,t) would 
vanish, and so its transform should vanish out there as well (i.e., both T and 6 
should vanish as x — oo). Therefore, the coefficient A p must vanish identically. 
Thus, 


0{x,p) = B p e ^ . (A-ll) 

Taking the transform of Eq. (A-7), one finds 

fdB\ 

UJ„ 0 = ~ Q(P) • < A - 12 > 

wherein 


Q(p) = £ P 


K 


Use of (A-ll) in (A-12) yields the 5 p -coefficient as 


(A-13) 


Bp = 


Qip) 
Vp ‘ 


(A-14) 


Therefore, from Eq. (A-ll), the transform of the surface temperature (at x = 0) 
is 
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(A-15) 



In a well-defined problem, this expression should comprise a well-defined analytic 
function for Re(p) > 0. The time-dependent surface temperature itself is then 
obtained from the inverse Laplace transform as 

<a - i6) 

wherein the integration contour. C, is comprised of a line parallel to the imag- 
inary axis in the complex p-plane. along which Re(p) > 0. The contour C is 
taken to lie to the right of all singularities of (Q(p)/ y/p), and we note that the 
presence of y/p places a branch cut along the negative- p axis. (Recall that Q(p) 
is essentially the transform of the incident energy flux.) 

It also is very important to note that the integrand in ( A-16) is usually a well- 
defined analytic function of p when Re(p ) > 0, because the Laplace transform 
integrals are themselves well-defined and convergent when Re(p) > 0 with Re(p) 
sufficiently large. Consequently, for r < 0, one may close the C contour at 
infinity in the right half-plane ({Re(p) > 0) so that |exp(pr)| — ► 0 at large \p\) 
to find that, since there are no enclosed singularities, 

T(0. r) = 0 for r < 0 . (A-17) 

The integrand, considered as a function of the complex variable, p, is to be 
“analytically continued'- into the left half-plane. Now, after the initial time, for 
t > 0, one can only close the contour at infinity in the left half-plane (wherein 
|exp(pr)| — 0 at large |p|), thereby enclosing the singularities of the integrand. 
Thus for r > 0, one generally finds that T( 0, r) / 0. 

To obtain a definite result, one needs to know Q(p), and it is therefore 
convenient to utilize a specific and simple model for the time-dependence of 
the energy flux on the surface. One might, at first, think that this must be 
done because the “backwards” unfolding problem of inferring a time-dependent 
incident energy flux from a known time-dependent surface temperature response 
might not have a unique answer. That is, because the instantaneous surface 
temperature depends upon the entire past history of the incident energy flux 


33 



in the form (obtained from (A-16) and (A-13) by means of certain Laplace 
transform relations), 


r(0,r)=^= f 

«\/x Jo 


q(T')dT' 

yjT — t' 


(A-18) 


it might be that several distinct irradiation histories could produce the same 
temperature response. (In more detail, we replace Q(p ) in (A-16) by the explicit 
transform, (A-13), interchange the order of integrations and note that p -1 / 2 
and (nr ) -1 / 2 constitute a transform pair.) This past-history integral represents 
the surface temperature response of a half-space irradiated at the surface by a 
time-dependent heat flux, q(r'). However, the inversion of Eq. (A-18) proves to 
be well-defined and unique! One finds (see. for example, The Mathematics of 
Physics and Chemistry, by H. Margenau and G. M. Murphy. D. Van Nostrand 
Co. (1956); Chap. 14), 


q(r 


_ jla. r 

yfic dr J 0 


Tj 0. r'lrfr' 
\fr — P 


(A-19) 


wherein we xecall that r = .D t h t. Thus, one could, in principle, infer the time- 
dependent energy flux incident on the surface from sufficiently accurate surface 
temperature data. 


The existence of a unique inversion to Eq. (A-18) is a nontrivial result, for it 
implies that there exists no time-dependent function (except "zero”), say Aq(r'), 
for which 


r A g (r')dr' _ Q 

Jo 

for all positive values of t. 

In the present experimental situation, though, the temperature data proba- 
bly are not known well enough to be practically utilized in a singular integration 
followed by a differentiation, as in Eq. (A-19). Therefore, we shall revert to 
solving the "forward” problem by assuming a reasonable, parameterized model 
for the time-dependence of the incident energy flux, and thereby shall calculate 
the surface temperature response. The parameters of the incident energy flux 
then are to be adjusted so that the computed temperature response agrees with 
the observations. 
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Our algorithm will be based upon two assumptions. First, the time depen- 
dence of the incident energy flux will be assumed as triangular, with a linear rise 
time t r to a peak flux, q max , immediately followed by a linear decay over a time 
duration t f . Second, the “rise” time. t r , and the "fall” time, tf, are to be taken 
from the time dependence of the power pulse into the CTX plasma gun, which 
itself proves to be roughly triangular. Thus, the only adjustable parameter in 
this simple but convenient model is qm*x- Accordingly, 

q(r) = for (0 < r < r r ) (A-20-a) 

q{r) = <?max (l - -f^) for (0 < r - t t < 77) (A-20-b) 

q(r) =0 for (77 < r - r r ) . (A-20-c) 

The Laplace transform of q{t)/n as defined by Eq. (20), may be calculated 
by utilizing simple integrations and elementary algebra. The result is 


e -p(r r +T } ) | 

P 2r f J ' 
(A-21) 

This Q{p) is a well-defined analytic function for Re(p) > 0. as is also Q(p)/^/p. 
The singularity of Q(p)/y/p is at the origin of the complex p-plane, and there 
is a branch cut along the negative-p axis. The inverse transform that yields the 
time-dependent temperature excursion on the tungsten surface is accordingly 
given by Eq. (A-16). Now, Eq. (A-21) contains three terms. Each of the three 
resulting contour integrals from (A-16) evidently has the form 

1 f e pT i 

Integral Term(j) = — / -rjrdp , j = 1, 2, 3 , (A-22) 

2tti Jc p 

in which the contour of integration, C, is parallel to the Im(p) axis and lies 
just to the right of the origin (Re(p) > 0 and Re(p ) — 0). Here, the t : ' s are 
obviously defined by (see Eq. (A-16)) 

"1 = r : t 2 = T - t t ; t 3 = T - T r - Tf . (A-23) 

Now, in all three cases, these integrals respectively vanish identically unless 
their respective Tf s are positive. Tj > 0. in view of the discussion just preceding 
Eq. (A-l 7). For times such that Tj > 0. the contour C may be deformed to lie 


q(p)= r 

Jo 


^Il e -Pr dr = ^ 




P 1 2 T t 


1 1 

1 

T r ^ 
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just below the negative-p axis, to pass around the origin in a counter-clockwise 
manner, and thence to pass again to infinity lying just above the negative-p 
axis. Making the change of variable of integration from p to p' = pr , . Eq. (A-22) 
becomes 


Integral Term(j) = rf jf jp^ip' ■ (A-24) 

Finally, changing the variable of integration from p ' to p" = (— p / ), and cor- 
respondingly changing to the appropriate integration contour in the p"-plane, 
the factor multiplying rf 2 may be trivially identified to be 1 /T(5/2). (See, for 
example, Modern Analysis, by E. T. Whittaker and G. N. Watson, Cambridge 
University Press (1952); Chap. 12, esp. p. 245.) Therefore, the use of (A-21) in 
(A-16) provides the following expression for the surface temperature excursion. 


T(0,r) 



(A-25) 


Here, u(t) is the unit step function (Heaviside function) that vanishes when its 
argument is negative and is unity when its argument is positive. The purpose of 
these step functions is to ensure that none of each of the three terms contributes 
anything unless its own Tj is positive. 

Recall that the time-like variables in (A-25) have been defined as t = D^t, 
where t is the actual time variable, and D t h = n/pC. At this point, it is 
convenient to redefine the time-like variables so as to be normalized to the rise- 
time t r . Thus, we shall define a new t , and a new Tj by 


T 3 tjt T ; Tj = tj/t r . (A-26) 

In terms of these new time-like quantities, Eq. (A-25) may be re-expressed as 
follows. 


n o,r) 


9max 

Wf) 



T 3/2 -u(t - 1)(1 + T f 1 )(t- 1 ft 2 

+ u(t - 1 - Tj)Tj l (T - 1 - ry) 3/,2 | 


(A-27) 
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Finally, when one notes that T(|) = Eq. (A-27) becomes the same as 

the expression for the surface temperature excursion that is given in the main 
body of this report as Eq. (4.1), together with Eqs.(4.2) and (4.3). 

The maximum temperature excursion thus occurs at the time when the fac- 
tor, {...} in Eq. (A-27), is maximized. Denoting {...} by /(r, 77), it is easy 
to see that /(r, 7y ) has a single maximum, which occurs when 1 < r < 1 4- Tf. 
The specific time of occurence of this maximum can be found as r = r P ^ = 
(1 + Tf) 2 /[(1 + Tf) 2 — 7y], which is plotted in Fig. 9 as a function of Tf. The value 
of /(r pea k, Tf), that is, the maximum value, is plotted in Fig. 8 as a function of 



Appendix B 
Plasma Heating 


B-l: INTRODUCTION 

In this Appendix we shall sketch several processes that may contribute to 
plasma heating in the CTX device. A recurrent theme shall be that the ions re- 
ceive the primary heating in the very process of ion production and entrainment, 
which they then share with the electrons by classical equipartition. However, it 
is recognized that certain supplementary processes involving microinstabilities 
may also play significant roles in electron heating, and these will be discussed 
in a later write-up. The consequences of these hypothetical descriptions prove 
to be consistent with experimentally inferred energy fluxes to the anode, with 
measured electron temperatures, and with measured neutral particle transverse 
energies in the plume. 

B-2: ION ENTRAINMENT HEATING 

Let us now consider the process of ion heating associated with ion creation 
and entrainment. To set the stage, we shall focus on some simple examples. 

First, consider two identical blocks, an “active” block A of mass M and 
a “passive” block B of mass M, each free to move on a common linear track 
without friction. The active block A contains a perfect motor that guarantees 
the constant velocity of block A in all circumstances, say va in the laboratory 
frame. (Alternatively, one could imagine that only A is rigidly attached to a 
constant-velocity conveyor belt.) Suppose block B is initially at rest (in the 
laboratory frame) in the path of oncoming block A. As viewed from A, B is 
approaching with velocity magnitude va- Let us now consider two kinds of 
collision of these two blocks, a perfectly elastic collision and a perfectly inelastic 
collision. 

In a perfectly elastic collision, all of the elastic potential energy stored in the 
compressed material during the collision of A and B is returned to the kinetic 
energy of the center of mass of B , which then recedes from the (motionless) block 
A with velocity va • In the lab frame, A continues with the constant velocity 
v A , whereas B now moves ahead of A with velocity 2v A . In the lab frame the 
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total momentum of the two-block system has tripled, and the total energy of 
the system has quintupled (thanks to the work done by the motor during the 
collision). 

However, in a perfectly inelastic collision, block B acts like a blob of putty, 
which sticks to A after the collision. Considered in the frame of A, all of the 
elastic potential energy that should have been built up during the collision, 
\Mv\, becomes distributed instead into many internal degrees of freedom, i.e., 
heat energy, and therefore none is available to be returned to the center-of-mass 
motion of B relative to A . In the lab frame, the momentum of the system has 
doubled, the kinetic energy of gross motion has doubled, and an amount of heat 
energy has been created that equals the additional kinetic energy. 

This simple example of a two-body collision illustrates the point that when 
matter at rest is entrained into a constant-velocity system, it must be heated by 
an amount equal to the additional kinetic energy that is created. 

Another example of this process is well known from ordinary gas dynamics. 
In the process of a strong planar shock impinging upon a compressible gas wdiich 
is at rest in the lab frame, the conservation laws of mass, momentum, and energy 
across the shock can be utilized to show that the thermal energy density behind 
the shock, P/( 7- 1), equals the kinetic energy density behind the shock, pV 2 j 2. 
Here, P and p are respectively the post-shock gas pressure and mass density, and 
V is the post-shock gas velocity in the lab frame. (This result can be obtained, 
for example, from the high-Mach-number limit of the expressions displayed in 
Fluid Mechanics (2nd edition), by L. D. Landau and E. M. Lifshitz, Pergamon 
Press (1987), Chap. IX — Shock Waves.) 

The presence of a sharp interface is not essential for this equipartition of flow 
energy and thermal energy. For, consider a uniform gas which is being legislated 
(by a distributed motor) to maintain a constant and uniform velocity V, and to 
which mass is being added uniformly at a rate p per unit volume. The rate (per 
unit volume) at which momentum is being created is 


pv = f , (B-l) 

wherein “/” is the externally legislated force-density analogous to the “motor” 
of the two-block collision example. The rate (per unit volume) at which flow 
kinetic energy is being created is 
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(B-2) 


l -pV 2 = W K . 

But the rate (per unit volume) at which the motor is performing work is 

fV = pV 2 = W M = 2W K ■ (B-3) 

Thus, only half of the work done by the motor becomes converted into flow 
energy. The other half of the work done must therefore be converted into heat 
energy, since this is the only other energy channel that exists. 

Let us finally consider the process of ion heating from the viewpoint de- 
veloped above, especially within the context of the Los Alamos CTX device. 
This device, like all MPD thrusters, with or without applied magnetic fields, 
is basically a crossed-field plasma accelerator. The voltage applied across the 
coaxial electrodes, in combination with the extant self and applied magnetic 
fields, forces the ion guiding centers to move with the local (E x B/B 2 ) drift 
velocity. (Here, E and B are respectively the local electric and magnetic field 
vectors.) In the case of self-field thrusters, the ( E t /Bq ) guiding center drift is 
downstream; whereas, for thrusters with strong applied magnetic fields, it is 
primarily azimuthal ( E r / B z ) pending its conversion to downstream motion by a 
diverging nozzle. But, in either case, the actual ion trajectories themselves are 
well known to be basically cycloidal, that is, having a circular gyration at the 
ion gyro-frequency superposed about the motion of the guiding center. (The ion 
is (cB/to,) in MKS units, for an ion of charge e and mass ui,.) 

Let us now consider the special character of such cycloidal orbits belonging to 
ions that have been created from resting neutrals (or from neutrals with velocities 
small compared to ( Ej B )). Such ions find themselves to be suddenly entrained in 
a moving electrodynamical system characterized by the (E/B) velocity. The only 
cycloidal orbits that satisfy the initial condition that the ions must start from 
rest are those that resemble a point on the rim of a coin standing on edge rolling 
on a table. Such orbits have the trivially obvious but very important property 
that their gyro- velocities are equal in magnitude to their common guiding center 
drift velocities. For a local ensemble of ions created at random from resting 
neutrals, these gyro-orbits have random phases and therefore can be regarded 
as comprised of random thermal velocities superposed on the flow velocity. In 
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the absence of frequent collisions, one has two degrees of freedom, hence 7 = 2 . 
The thermal energy density of the ions is thus 

(7 - l) -1 ^ = n,Ti = y-n, = ~niv 2 thl , (B-4) 

which is the same as the flow kinetic energy density for these orbits since v t h, = 
E/B. (Here P t is ion pressure perpendicular to B , n is ion number density, 
T{ is ion temperature, and Vthi is ion “thermal velocity.”) Although one has a 
cylindrical-shell distribution function rather than a Maxwellian, the “thermal” 
pressure effects across B are close enough that the analogy is useful, and Eq. (B- 
4) can be regarded as defining the ion temperature for these orbits wherein v t hi = 
E/B. In fact, the exact Eq. (B-4) can be directly derived for the ion pressure 
associated with a cylindrical shell distribution function /;( u) = (ni/ 7 r)S(u± - 
v thi)6( u z), where the velocity vector u represents the difference between the 
particular velocity and the average velocity, with the local z-direction along B. 
Thus, in this generalized sense, the thermal energy density of ions, created from 
resting neutrals in a crossed-field plasma accelerator, equals their flow energy 
density. (Note that it was not necessary to invoke collisionality here, which serves 
to emphasize just how fundamental is this equipartition between flow energy and 
thermal energy during entrainment heating.) 

All of the above examples illustrate the fundamental process of entrainment 
heating, but to complete this discussion of ion heating in the CTX device we 
must now deal with those thus-far-suppressed topics of electron-ion interactions, 
ion viscous heating, and ion thermal conduction losses. These processes, in con- 
juction with the inelastic process of ion entrainment, all contribute to setting the 
level of the ion temperature. These topics will now be dealt with in connection 
with 40 MW class operation of the CTX device. 

By taking moments of the Boltzmann equation for the ion phase-space distri- 
bution function, while allowing for the creation of ions from resting neutrals, one 
can obtain the following equation governing the evolution of ion temperature, 
Ti. 


1 dT, 
7 — 1 eft 


+ T,V ■ Vi 




+ ••• 


(B-5) 


Here, Vi is the ion flow velocity, ( h/n ) is the local ion production rate per ion 
(usually dependent on the electron temperature, T e ), ^ is the total convective 
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derivative, (jj + Vi ■ v), and the dots (■ • •) signify the additional complicating 
effects of electron-ion equipartition of thermal energy, ion viscous heating, and 
ion thermal conduction losses. We shall presently discuss each of these topics as 
it influences the ion temperature. (The thermal source term, m i V i 2 /2, eventually 
disappears downstream by the time neutrals may have been accelerated up to 
the ion flow velocity.) 

But first, it is of interest to remark that if none of these complicating effects 
are dominant, then the mathematical structure of Eq. (B-5) already provides 
to us information that supplements the simple guiding center model described 
earlier. Namely, if V • also does not dominate, then Eq. (B-5) indicates the 
tendency for the ion temperature to "settle in” at the level of the ion flow energy 
(per particle), in accordance with the simple model of entrainment heating; and 
moreover that the relaxation time for a deviated T,- to approach this special 
value is determined by the ion production rate (per ion). (See also H. A. Hassan, 
Physics of Fluids 9 (1966) 2077.) Thus, Eq. (B-5) essentially verifies the inelastic 
heating and the equipartition of ion flow energy and ion thermal energy due 
to the ion production and entrainment process, and it yields a characteristic 
time scale for such heating; but it also confronts us with modifications and 
qualifications to the simple picture that have to be considered. 

With regard to electron-ion temperature equipartition in the CTX device, 
it can be observed here that, if the ion temperature is indeed near the flow 
energy, then calculations to be presented in a subsequent write-up show that 
equipartition with electrons is, in fact, a modest effect, thus being a posteri- 
ori consistent with the provisional neglect of this effect in Eq. (B-5). But even 
in the presence of rapid temperature equipartition between ions and electrons, 
conclusions can be drawn similar to those already presented here. An equation 
assuming rapid electron-ion equipartition (but neglecting the Hall effect and as- 
suming that the plasma is optically thick to excited-state decay radiation) has 
been derived elsewhere that is a generalization of Eq. (B-5) to include excitation 
and ionization processes in hydrogen. (See ^Preliminary Scoping Studies for 
Nozzle-based Coaxial Plasma Thrusters,” by R. A. Gerwin, K. F. Schoenberg, 
and D. J. Rej; in Eighth Symposium on Space Nuclear Power Systems. Albu- 
querque, 1991; Proceedings Pub. by American Institue of Physics, AIP Conf. 
Proc. 910116; Part Two, p. 493.) This equation can be expressed schematically 
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as 


2 dT - / 1 , 2 \ 

— + 2TV * Vj = (effective rate) ( -m t - V 2 -T — A£ H . (B-6) 

7 - 1 at \2 7 — 1 / 

wherein T is the common electron-ion temperature and (• • •) again refers to vis- 
cous and resistive heating and thermal conduction loss terms. (Near to a steady 
state, it can be demonstrated that the effective production rate is practically 
the rate of electron impact generation of excited hydrogen states.) The impor- 
tant new element in Eq. (B-6) is the appearance of A E, an energy increment 
characteristic of electron-impact-induced atomic transitions that also generally 
depends on the relative population levels and T-dependent rate coefficients. In 
hydrogen near to a steady state, it can be shown that A E will be of order 10 
eV. The mathematical structure of Eq. (B-6) contains the features of a “critical 
ionization velocity” process, in the sense that if the flow energy is too small, 
T will fall and not produce more ionization; whereas if the flow energy is suffi- 
ciently large, then T can rise to a level that continually produces ionization. (See 
H. Alfven, “Collision Between a Non-ionized Gas and a Magnetized Plasma,” 
Rev . Mod. Phys. 32 (1960) 710.) In fact, if the flow energy is sufficiently large 
against A £, then Eq. (B-6) implies, as before, that “the thermal energy will sat- 
urate near the flow energy.” In the present case of rapid equipartition, however, 
it is the thermal energy of an electron-ion pair that would saturate near the flow 
energy of an ion. 

In the CTX device at 40 MW operation, guiding-center drift energies in the 
breech are inferred to be about 25 eV, and axial neutral energies in the plume 
have been measured at upwards of 100 eV. In view of the above discussion, 
it may therefore not be unrealistic to expect “ion entrainment heating” at a 
thermal level of at least 25 eV in these experiments. However, the production 
and sustainment of such an ion temperature is subject to consideration of ion 
viscous heating and thermal conduction losses, upon which we now embark. 

B-3: ION VISCOUS HEATING AND THERMAL 
CONDUCTION COOLING 

Within the context of 40-50 MW operation of the CTX device, measured 
densities at the muzzle(~ 1.0 x 10 2 ° m~ 3 ), self-magnetic field strengths (6 — 
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8 x 10 2 tesla), and 25 eV or higher ion temperatures according to the ion 
entrainment model, when taken all together imply that the deuterium ions are 
well-magnetized. 


Ueim ~ 10 . 

(Here, «*>„• is the ion gyro-frequency and r tt is the ion-ion collision time. The 
ion-ion collision frequency is written as u {i - T~ t l .) Moreover, since the axial 
transit time is only a few ^ s (based upon the measured axial velocity of neutrals 
of about 10 5 m/s and upon the indication of electrode wear along a 10-30 cm 
axial distance), there can be only a few ion-ion collisions during an axial transit 
time t z . 


t z va ~ a few. 

Therefore, a strongly-magnetized, almost collisionless model for the ion behavior 
is appropriate, provided also that ion-neutral collisions are negligible. (There 
is, at present, no information regarding the density of neutrals during the CTX 
discharge.) Such a model for the ions shall be used here to estimate the effects of 
ion viscosity and ion thermal conduction. A very detailed and complete reference 
on this subject, transport coefficients of fully ionized plasmas, exists. (S. I. 
Braginskii, Transport Processes in a Plasma.'’ in Reviews of Plasma Physics , 
Vol. 1, Consultants Bureau, M. A. Leontovich, Ed., New York (1965).) We shall 
use the results of Braginskii for strongly magnetized ions. 

Suppose there exists a radial gradient (width S r ) of axial velocity (V z ), due to 
a no-slip boundary condition. Then, one can estimate that the fractional increase 
of ion temperature due to viscous heating during an axial transit (length A*) is 
given by 

°- 3 (U^ (x 1 ) (x) ( viscous ion heating). (B-7) 

For {udTii) 2 ~ 100, V 3 ~ 10 7 cm/s, A* ~ 10 cm (nominal), S r ~ 2 cm ( 6 r ~ half 
of the channel half-width and also two gyro-radii), tu ~ 2 x 10 -6 s. the above 
expression becomes 


(viscous ion heating) ~ 0.15 . 


(B-7-a) 



Suppose there exists a radial gradient (width <5 r ) of ion temperature T t . Then, 
one can estimate that the fractional loss of ion temperature due to cross-field 
thermal conduction during an axial transit is given by 

ru ( 2T t /m t )A z 1 ( vtkiTji \ f vthi A ( 

(UciTii) 2 S 2 r V 2 "(U ^) 2 V S r J \ V z J \ S r J 

(thermal conduction ion cooling). (B-8) 

Using the same parameters as above, and with v t hi = 0.5 X 10 ' cm/s correspond- 
ing to E/B in the breech, this expression reduces to 

(thermal conduction ion cooling) ~ 0.13 . (B-8-a) 

The first thing to notice about these estimates is that they are each small frac- 
tions of the estimated ion temperature; secondly, they are almost equal in mag- 
nitude; and thirdly, these two processes tend to cancel one another. We therefore 
conclude that the simultaneous omission of these two processes in the consider- 
ation of Eq. (B-5) is justified. 

Now, these two numerical estimates depended upon assumptions in CTX 
regarding gradient scale lengths, and axial transit distances; and they w’ould be- 
come individually larger in regions of sharper gradients (smaller 6 r 's) or larger ax- 
ial transits (larger A^’s). However, both expressions scale identically, as A 2 6~ 2 , 
so that their mutual cancellation would be retained. Moreover, 6 r itself cannot 
be much smaller than already assumed, since the ion gyro-radius under the given 
conditions is about 1 cm. 

Note that the ratio of viscous heating (B-7) to conduction cooling (B-8) 
is about Vf /(3v 2 hi ). Since V z represents a flow velocity in the nozzle ( Va) 
which has begun to be accelerated beyond the drift velocity in the breech, one 
might expect that this ratio is always of order unity — as can be experimentally 
inferred. (We are tacitly assuming here that v t ki represents the ion temperature 
of particles entrained upstream that have been swept into the nozzle region from 
the breech. This assumption is a posteriori consistent with the estimate that 
net heating or cooling effects on the ions during their journey from breech to 
nozzle are small-to-moderate.) 
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Appendix C 

Inframetrics IR Camera Calibration 


The Inframetrics model 600 IR video camera measures temperature change on a 
target surface by measuring changes in infrared emission from that surface over 
a 9 to 11 micron bandwidth. The infrared intensity reaching the camera is the 
sum of radiation emitted by the target and the background radiation reflected 
off the target. Here, the radiation emitted by the target is the Planck black-body 
function integrated over the bandwith and weighted by the target emissivity e. 
The reflected radiation is related to the background temperature weighted by 
the target reflectivity (l-£). 

In many applications, the background temperature is uniform and static in time. 
Hence, an absolute conversion of the emitted IR radiation into temperature sim- 
ply requires a value for the target emissivity and the background temperature. 
During CTX discharges, however, all electrode surfaces are heated by the plasma. 
Hence, the nominal background temperature is both time dependent and near 
that of the imaged electrode surface. For this case, the camera “emissivity pa- 
rameter” and “background temperature parameter” are ill-defined. 

Absolute calibration of the IR camera, for CTX relevant conditions, entailed 
uniformily heating the thruster to 65° C using hot air. The camera was then 
positioned to replicate the spatial line-scans used to acquire the reported data. 
As the electrodes slowly cooled, covering the range of data-relevant temper- 
ature changes, true electrode temperature was monitored by a thermocouple 
attached to the anode (at the position of the line-scan). The camera “emissiv- 
ity parameter was then continuously adjusted to force agreement between the 
camera inferred temperature and that of the thermocouple. Here, the camera 
“background temperature parameter” was set at 20° C. In this way, a calibra- 
tion “curve” was generated that permitted conversion of imaged IR intensity 
into electrode temperature. This calibration procedure did not account for the 
effects of differences in electrode temperature or spatially distributed electrode 
temperature. Both these effects are considered small for the CTX data reported. 
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Table I. 

Parameters and Peak Energy Flux Estimates for 40-MW-Class Shots 

Shot # 

U(ps) 

Peak Power 
(MW) 

tj(ps) 

(10% of 
Peak Power) 

T f 

fm ax 

(l + -13r/) 

(AT) m ax 

(°C) 

Qmax 

20274 

100 

50 

330 

3.30 

1.43 

11.7 

23 

20280 

100 

60 

300 

3.00 

1.39 

7.2 

14 

20291 

100 

40 

420 

4.20 

1.55 

13.5 

24 

20300 

90 

45 

400 

4.44 

1.58 

11.7 

22 

20301 

100 

50 

340 

3.40 

1.44 

15.3 

30 

20338 

100 

40 

450 

4.50 

1.59 

12.6 

22 

20340 

75 

48 

420 

5.60 

1.73 

12.6 

23 

20341 

110 

50 

360 

3.27 

1.43 

8.1 

15 

20342 

110 

50 

360 

3.27 

1.43 

9.9 

19 

20373 

100 

43 

425 

4.25 

1.55 

13.5 

24 

20375 

100 

46 

400 

4.00 

1.52 

15.3 

28 

20377 

85 

50 

440 

5.17 

1.67 

18.9 

34 


Assuming (AT) (°C) at t = 0.5 ms is (AT) max (°C). 
mean power: p = 48 ± 5 MW 
mean: q max = 23 ± 5 [MW /m 2 ] 
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Table IA. Peak Energy Flux Estimates for 40-MW-Class Shots 


Shot # 

r(500/rs) Peak Power (MW) 

T f 

/(t "/) 

W/m 2 ) 

20274 

5.0 

50 

3.30 

0.93 

35 

20280 

5.0 

60 

3.00 

0.85 

24 

20291 

5.0 

40 

4.20 

1.28 

30 

20300 

5.6 

45 

4.44 

1.18 

29 

20301 

5.0 

50 

3.40 

0.96 

45 

20338 

5.0 

40 

4.50 

1.40 

25 

20340 

6.7 

48 

5.60 

1.31 

31 

20341 

4.5 

50 

3.27 

1.03 

21 

20342 

4.5 

50 

3.27 

1.03 

26 

20373 

5.0 

43 

4.25 

1.30 

29 

20375 

5.0 

46 

4.00 

1.18 

36 

20377 

5.9 

50 

5.17 

1.39 

41 

Assuming (AT) (°C) at t = 

0.5 ms is not (AT) mai 

(°C). 




2.8 AT[°C] 
-(lOO^Xs] f(T,Tf) ' 


with t = 


500ps 

t r 


and AT (°C) measured at 500 p s. 


mean power: p = 48 ± 5 MW 
mean: q max = 31 ± 7 [MW/m 2 ] 


Table II. Parameters and Peak Energy Flux Estimates for 10-MW-Class Shots 

Shot # 

t T (v s) 

Peak Power 
(MW) 


T S 

fmax 1 

(AT) mai 


20381 

110 

12 

400 

3064 

1.47 

3.6 

6 

20382 

120 

9 

320 

2.67 

1.35 

9.0 

17 

20383 

100 

15 

250 

2.50 

1.33 

7.2 

15 

20384 

100 

9 

380 

3.80 

1.49 

8.1 

15 

20385 

100 

10 

350 

3.50 

1.46 

6.3 

12 

20386 

100 

10 

360 

3.60 

1.47 

6.3 

12 

20387 

100 

11 

330 

3.30 

1.43 

4.5 

9 

20388 

110 

11 

350 

3.18 

1.41 

6.3 

12 

Assuming (AT) (°C) at t = 0.5 ms is 

(A T) max (°C). 



mean power: p = 

11 ± 2 MW 






mean: q max = 12 ± 3 [MW/m 2 

] 
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Table IIA. Peak Energy Flux Estimates for 10-MW-Class Shots 
((AT) max not assumed to be at O.o ms) 


Shot # 

t( 500/xs) Peak Power (MW) 

T f 

f(T Tf) 

<?ma*(MW/rn 2 ) 

20381 

4.5 

12 

3.64 

1.20 

8.0 

20382 

4.2 

9 

2.67 

0.88 

26 

20383 

5.0 

15 

2.50 

0.72 

28 

20384 

5.0 

9 

3.80 

1.10 

21 

20385 

5.0 

10 

3.50 

1.00 

17 

20386 

5.0 

10 

3.60 

1.03 

17 

20387 

5.0 

11 

3.30 

0.93 

14 

20388 

4.5 

11 

3.18 

1.00 

17 

Assuming (AT) (°C) at t = 

0.5 ms is not (A T) max 

(°C). 



500/is 


2.8 Ar[°C] 

/t r [100/xs] f{r,T f ) 

and AT (°C) measured at 500 ps. 


mean power: p = 11 ± 2 MW 
mean: q max = 18 ± 5 [MW/m 2 ] 
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SCALE - cm 


(a) Schematic. 


(b) Photograph of the CTX. 


(c) Photograph of the CTX facility. 


Figure 1 .—Coaxial plasma thruster. 


ORIGINAL PAX. 

SLACK AND WHITE PHOTOGRAPH 
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microseconds 

Figure 2. — Nominal "round-top" voltage and current waveforms 
for a 40 MW discharge. 


Figure 4.— Infrared image of the CTX coaxial plasma thruster. 


Position (m) 

Figure 3. — Vacuum field calculation of applied magnetic 
field structure at the time of discharge. 




Normalized Surface Temperature Rise f(x, t ( ) 


0 



(a) Space-time resolved surface temperature rise for 
discharge 20383. 



(b) Position of IR spatial scan line for discharge 20383. 
Figure 5.— Discharge 20383. 



Normalized Time 


Figure 6. — Time-dependence of normalized surface temper- 
ature rise for T fa j| = 5.5 T hse . 



Figure 7.—' Time-dependence of normalized surface temper- 
ature rise for T faU = 3.0 T rise . 



Microseconds 


Figure 8. — Data showing surface temperature vs. time. Each 
narrow vertical strip is a compressed spatial line scan. 



Figure 9.— Normalized peak surface temperature rise vs. 
( T fal/ T rise) ~ T f* 
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Plasma Potential (Volts) 



Figure 11 . — Plasma potential vs time: (a) probe in muzzle, 1 mm from anode, (b) probe in muzzle, 5 cm from anode, (c) probe 
40 cm downstream of thruster exit 
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Spectral Density 








total Plasma Potential (V 



Figure 14. — Plasma potential within the muzzle (also shown Is the plasma potential signal following a 20 kHz low pass filter) 
and fraction of total discharge current flowing through the tank wall. 



0 200 400 600 800 1000 

Mteroseconds 

Figure 1 5. — Input power vs. time for 1 0 MW discharge 20363. 
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Figure 1 6. — Radiation profile shows twin peaks across the anode-cathode-anode gaps. 



Time, ^s 


Figure 17. — An arrival time distribution from the TOF neutral 
particle spectrometer for 40 MW magnetized coaxial gun 
operation. The temporal difference between the two signal 
peaks represents the transit time of deuterium neutrals dawn 
an approximate 2 meter drift tube. 
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